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Abstract An overview of the vector and tensor Generalized Parton Distributions for 
a charged pion is presented. Such observables, belonging to the set of quantities funda- 
mental for a detailed investigation of the hadronic inner dynamics, have been evaluated 
within a fully covariant Constituent Quark Model, based on a proper Ansatz of the 
pion Bethe-Salpeter Amplitude and the Mandelstam formula for matrix elements of 
operators acting on relativistic composite systems. Given the very encouraging results 
already obtained for the vector distribution of a charged pion, the model has been 
extended to the tensor one, and some preliminary calculations will be illustrated. 

Keywords Pion Generalized Parton Distributions • Covariant Constituent Quark 
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1 Introduction 

The non perturbative regime of Quantum Chromodynamics (QCD), namely the present 
theory of the strong interaction, is still a paramount challenge. In particular, achieving 
a complete 3D description of hadrons is the topic of many and coherent efforts, both 
on experimental and theoretical sides. The perturbative regime, with its fundamental 
feature, the asymptotic freedom, has been experimentally investigated in great detail, 
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and this has allowed to access the short-distance behavior of the hadronic wave func- 
tions. But, in order to shed light on confinement, the other peculiar feature of QCD 
(deeply related to its non linearity), one has to move towards the external part of the 
hadron wave function. 

From the theoretical point of view, the Light-front (LF) framework, with variables 
defined by: a = a ± a 3 and aj^ = {a x , a y }, is very suitable for describing relativistic, 
interacting systems, like hadrons. Among the several motivations for adopting the LF 
framework, beyond the well-known fact that the dynamics of the light-cone is naturally 
described in terms of LF variables, one has to mention that the vacuum is almost trivial 
within a LF field theory [T] . Therefore, e.g. for the pion, one can construct the following 
meaningful Fock expansion (with some caveat for the zero-mode contributions [T]) 

k) = \qq) + 1 95 95) + 1 99 g) 

where the valence component \qq) and the non valence ones are fully considered. Within 
our approach, based on a covariant description of the pion state (see [2 , 3 and references 
quoted therein), one can explore also the contributions beyond the valence term. 

In recent years, it has been recognized that a wealth of information on the partonic 
structure of hadrons is encoded in the Generalized Parton Distributions (GPD) (see, 
e.g., Ref. [1]), as well as in the Transverse-momentum Distributions (TMD) (see, e.g., 
Ref. [5]). Those observables, that can be experimentally accessed through Deep Virtual 
Compton Scattering and Semi-inclusive Deep Inelastic Scattering, respectively, are 
expressed by quantities invariant for LF-boosts, and allow one to parametrize matrix 
elements (between hadronic states) of relevant quark-quark or gluon-gluon correlators. 
Aim of our covariant, phenomenological analysis is the evaluation of the leading-twist 
pion GPD's. 

2 Generalized Parton Distributions 

The generalized parton distributions are off-diagonal (respect to the hadron four- 
momenta, i.e. pf p{) matrix elements of quark-quark (or gluon-gluon) correlator 
projected onto the Dirac basis (see, e.g., Ref. [6] for a thorough investigation of the 
pion case), and naturally summarize a lot of information, scattered in many observ- 
ables, like e.g. electromagnetic (em) form factors or parton-distribution functions. In 
particular, the pion, i.e. the simplest hadron to be considered, has two leading-twist 
quark GPD's: i) the vector, or no spin-flip, GPD, H^(x, £, t), and ii) the tensor, or spin- 
flip, GPD, j*(x, £, t) (I = IS, IV labels isoscalar and isovector GPD's). In order to 
avoid Wilson-line contributions, one can choose the light-cone gauge [4] and gets 
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^ , " )/, (p'l?,(^)^ U(~*)|P>| 5=0 (2) 

where z = {z + — z° + z 3 , z±}, 4>q(z) is the quark-field isodoublet and 

- = ^ . e= - ^ . ^-^ 2 - . ^ = ^ (3) 

with {fc + — Z\ + /2, kj^ — Aj_/2} the initial LF momentum of the active quark. Notice 
that the factor of two multiplying the vector GPD is chosen in order to obtain the em 
form factor from J dx H(x,£,t) (see below). 

As anticipated above, there exist basic relations between the GPD's and other 
relevant observables. If we consider a charged pion, the vector GPD is related both to 
the parton distributions, q(x), and to the em form as follows 

H%(x, 0, 0) = 6{x)u{x) - 9(-x)u(-x) 



F 7r (t)= J dx Hi v (x,t;,t)= J dxH%(x,t,t) (4) 

where H% — H^ S + . The relation between the non-spin flip GPD and the em form 
factor of a charged pion can be generalized, if one considers Mellin moments of both 
vector and tensor GPD's, obtaining the so-called Generalized Form Factors (GFF). In 
particular, for the it-quark vector and tensor GPD's, the corresponding Mellin moments 
are given by 

/ dxx n H%{x,i,t)=Y J (n) i Al +lj2i {t) , 
i 

dxx n E% >T (x,Z,t) = Y^^Bn+lMt) (5) 



where I — {(n + l)/2} ({a} is the integer part of a), A^ +1 2 i(i) is a vector GFF and 
•®n+i 2iW a tensor one. One has isoscalar (isovector) GFF's if n is odd (even). It is 
worth noting that 



l 

dxH%(x,£,t)=Al (t) = Fi T (t) (6) 

l 



and 

el 

dxE^ T {x,i,t) = Bla(t) (7) 

l 

with £>"q(0) 7^ 0, the tensor charge for n = 1, or the tensor anomalous magnetic 
moment (as it is called in Ref. [7])- 
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In order to attach a physical interpretation to the GFF, one can generalize to a 
relativistic framework what is familiar for the em form factors (scalars and then boost- 
invariant observables) in a non relativistic approach. In this case, one interprets the 
3D Fourier transforms of the em form factors of a hadron as the intrinsic (Galilean- 
invariant) em distributions in the coordinate space (e.g., for the pion, one has the charge 
distribution, while, for the nucleon, one has both charge and magnetic distributions). 
In the relativistic case, one has to consider Fourier transforms of GPD's, that are 
quantities depending upon variables invariant under LF boosts. Moreover, only the 
transverse part of A^ can be trivially conjugated to coordinate-space variables, while 
for x and £ (proportional to A + ) this is not possible. Therefore one can introduce 
bidimensional Fourier transforms with respect to Aj_, for gaining some insight in the 
spatial distributions of the quarks, still keeping the description invariant for proper 
boosts (i.e. LF boosts). In particular, from Eq. |5| one can see that, putting £ = 0, 
only A^ +1 o(^ 2 ) (°f a l so Eq. ^) and B" +1 q(A ) survive. There exist an infinite set 
of frame where £ = 0, due to the LF-invariance of £, that are indicated as the Drell-Yan 
frames. In these frames, where A + = and A 2 = 
transforms of the previously defined GFF's as follows 



Aj_, one can introduce the Fourier 



A n (b ± ) = 



W) 



(Ai-bj 



AnM A ) , B n (b ± ) = 



/ 



W) 



zA x b j 



B nfi (A z )(8) 



where bj^ = |bj_|. They have a direct physical interpretation as quark densities in the 
impact-parameter space (IPS) [HE]. In particular, A n {b±) represents the probability 
density of finding an unpolarized quark in the pion at a certain distance bj_ from the 
transverse center of momentum. If one takes into account the quark polarization, then 
the probability density of finding in the Drell-Yan frame a charged parton of fixed 
transverse polarization, sj_, and given bj_, is 



/On (b_L,Sj_) 



1 



where 



r n (b ± ) 



A n (b ± ) + r n (b ± ] 



1 d B n (b ± ) 



Im-r: d b^_ 



(9) 



(10) 



The IPS density p n can be obtained by calculating the proper matrix elements of the 
projector for a transversely-polarized quark (s^ = {0,0, sj_}), namely 



1 + H 1 1 

-7 [1 + s ± ■ 7_l7 5 ] = - 



7 +? 
S a 



(ii) 



It should be pointed out that the distribution density of a quark with a given 
helicity (longitudinal polarization) is expressed by the first term in Eq. (j9|, since the 
pion is a pseudoscalar meson and the quark helicity operator (for a massless fermion: 
75 /2) has a vanishing expectation value, due to parity invariance. 

Equation |9| clearly indicates how to access the quark distribution in the impact- 
parameter space through the study of the GPD's. Moreover, as a closing remark, one 
could exploit E^x to extract more elusive information on the quasi-particle nature 
of the constituent quarks, like their possible anomalous magnetic moments, once the 
vector current that governs the quark-photon coupling, is suitably improved (see below 
and [10] for a discussion in the lattice framework). 
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3 Transverse momentum distributions 



The transverse momentum distributions are diagonal (in the pion four-momentum) 
matrix elements of the quark-quark (or gluon-gluon) correlator with the proper Wilson- 
line contributions (see, e.g., Ref. [6]) and suitable Dirac structures. The TMD's depend 
upon x and the quark transverse momentum, kj^ (notice that = leads to £ = 
t = 0), and since the integration over kj_ is not performed, the Wilson-line effects 
must be carefully analyzed. At the leading-twist one has two TMD's, for the pion: 
the T-even fi(x, |kj_| 2 ), that yields the probability distribution to find an unpolarized 
quark with LF momentum {a;,kj_} in the pion, and the T-odd hj^(x, | k_j_ | 2 , 77) , related 
to a transversely-polarized quark and called Boer-Mulders distribution [TT]. The two 
TMD's allow one to parametrize the distribution of a quark with given LF momentum 
and transverse polarization, i.e. (see, e.g., Ref. [12] ) 

p(x,\t ± ,s ± ) = ^ 



/i(a;,|k ± i 2 ) + i 



hi(x, |k_L | 2 , 77) 



(12) 



where the dependence upon the variable r\ in hj~ can be traced back to the Wilson-line 
effects. Moreover, it is worth reminding that bj^ is not conjugated to k^, but to Aj_. 

By choosing the light-cone gauge and the advanced boundary condition for the 
gauge field, the effect of the Wilson lines (final state interaction effects) can be trans- 
lated into complex phases affecting the initial state (see, e.g., Ref. 13 ). At the lowest 
order, the unpolarized TMD // (x, |kjj 2 ), acquaints the following form 

|k x | 2 ) \ 



fl V (x,\k ± \ 2 ), 



1 



1 



^e^ + -/ 2 -^-](p|^(-^ )7 +rj^(^)|p>| z+=0 , (13) 
After integrating over k^, one gets the standard unpolarized parton distribution q(x), 



q(x) = jdk ± /«(x,|k x [ 2 ) = Hf(x,0,0) . (14) 

The T-odd TMD, h^(x, | k_i_ | 2 , 77) vanishes at the lowest order in perturbation theory, 
since it becomes proportional to the matrix elements 



<j# 9 (-2*)*<r™ I I M^)\p)\ z+=0 , (15) 



that are equal to zero, due to the time-reversal invariance. In order to get a non van- 
ishing Boer-Mulders distribution, one has to evaluate at least a first-order correction, 
involving Wilson lines (see, e.g., Refs. [14] and [6]). 
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4 The Covariant Constituent Quark Model 

Presently, extensive theoretical and experimental research programs are being pursued 
to gain information on both GPD's and TMD's of hadrons (particularly of the nucleon). 
Within such efforts, considering that the lattice calculations (in Euclidean space) can- 
not cover the whole kinematics that experiments can investigate, it seems worth con- 
sidering the development of phenomenological models, that are able to include general 
features, like the covariance, and to properly describe both the experimental data and 
the lattice ones. In this way one could have an effective tool for the analysis of the re- 
sults of the forthcoming experiments (e.g., the ones planned and approved at TJLAB), 
devoted to set-up the 3D tomography of hadrons. 

The main ingredients in our covariant constituent quark model (CCQM) are rep- 
resented by i) a model of the 4D quark-hadron vertex, namely the Bethe-Salpeter 
amplitude, and ii) the extension to the GPD's and TMD's of the 4D Mandelstam ex- 
pression [15], originally introduced for calculating matrix elements of the em current 
operator for a relativistic interacting system. It is important to emphasize that our 
investigation, naturally goes beyond a purely valence description of the pion [HE]- 



4.1 The Mandelstam Formula for the electromagnetic current 

The Mandelstam formula allows one to express the matrix elements of the em current 
of a composite bound system, within a field theoretic approach 15 . It has been applied 
for evaluating the form factors of both pion [TB] and nucleon [T7]. Furthermore, it has 
been exploited for calculating the vector GPD of the pion [2][3] and for a preliminary 
evaluation of the tensor GPD [18] ■ 

For instance, in the case of the em form factor of the pion, in the spacelike region, 
the Mandelstam formula within CCQM reads (see, e.g., Ref. [16]) 

X Tr[S(k - P„) j 5 S(k + q) r^{k, q) S(k) 7 5 ] (16) 

where f„ = 92.4 MeV is the pion decay constant, N c — 3 is the number of colors, 

m = 220 MeV the CQ mass, S(p) = the fermionic propagator, and P^ik, q) 

p — m + le 

the quark-photon vertex, that simplifies to 7^ in the spacelike region. In presence of 
a CQ, one could add to the bare vector current a term proportional to an anomalous 
magnetic moment, namely a term like 

2m q 

that should not be confused with the quark contribution to the nucleon anomalous 
magnetic moment. In the present approach we have disregarded such a possibility (cf 



|10j for an improved vector current within a lattice framework). In Eq. (16 1, it has 



been assumed a pion Bethe-Salpeter amplitude (BSA) with the following form 

W{k-P,p) = -— S(k-A/2) 7 5 A(k-P,p) S(k-P) (17) 
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where only the most important Dirac structure 75 has been retained (see Ref. [19] for 
a detailed discussion of terms beyond 7 s ), and A(k — P,p) is the momentum-dependent 
part. 

The generalization to the case of GPD's, can be found in Ref. [213] f° r the vector 
GPD, and in 18j 20] for the tensor one. In particular, both an accurate description 
of the em form factor of the charged pion for t < [2j[3] and a nice behavior of the 
parton distribution evolved up to the experimental scale (see 20 ) can be achieved by 
adopting the following analytic covariant Ansatz for the momentum-dependent part of 
the BSA 

A(k-P,p) = - ^ T (18) 

(k - A/2) 2 - m| + ie\ (P - fc) -m% + ie\ 

where the parameter mji = 1200 MeV is adjusted to fit f^, while the constants C is 
fixed through the form factor normalization, F^{t — 0) = 1, that acts as a normalization 
condition for the BSA. 

For the sake of concreteness, let us give both vector and tensor GPD's for the u 
quark, viz 

2H u (x,£,t) = -t1l j ^i s l p+x - k+ ] A{k-P,p') A{k-P,p) 

xTr js (k - P) 7 5 S (k + A/2) 7+ S (k - A/2) 7 5 | (19) 



and 



-Kt(x, £, t) = inj 44iS[P + x - k+] A(k - P,p') A{k - P,p) 



P+A 1 - P j A + „„ , ^ .„ f d 4 k 

P+mir 

xTr[S{k - P)j°S(k + z4/2) 7 + 7 J S(k - Z\/2) 7 °] (20) 



where 1Z = 2N c m 2 / f% and j = 1, 2. The 6 function imposes the correct support for the 
active quark, i.e. |£| < x < 1 (see 2,3 ), corresponding to the so-called DGLAP region 
|21j . The diagrammatic representation of the vector GPD, with a LF-time ordering, is 
given in Fig. la and lb. 



A A 




Fig. 1 LF time-ordered analysis of the pion GPD within the CCQM framework. Diagram 
(a): contribution of the active-quark in the valence region with 1 > x > Diagram (b): 
contribution from the pair-production mechanism in the nonvalence region, |£| > x > — 
(After Ref. [2151 1. 



X 



X 



Fig. 2 Preliminary isoscalar and isovector tensor GPD's for a charged pion, within CC QM , 
for positive x. The behavior for negative values of x can be obtained by recalling Eq. i |- 1 [ )■ 
Thick solid line: £ = and t = 0. Thick dotted line: £ = and t = -0.4 GeV 2 . Thick dashed 
line: £ = and t = -lGeV 2 . Thin dotted line: £ = 0.96 and t = -0.4 GeV 2 . Thin dashed 
line: £ = 0.96 and t = -1 GeV 2 . 



In the CCQM investigation of GPD's, a Breit frame with A + = — A~ > is chosen 
in order to access the whole range of the variable £, i.e. —1 < £ < 1. In this way both 
the valence and non valence region can be investigated, addressing the interesting topic 
of the smooth transition from the DGLAP (valence) regime to the ERBL (non valence) 
one [22] . In Ref. 0C5], the vector GPD and the unpolarized TMD, resulting from the 
CCQM approach, are thoroughly analyzed. In particular, the CCQM calculation was 
able to carefully reproduce the experimental data of the em form factor of the charged 
pion over the whole spacelike range of t, presently explored. As to the GFF's obtained 
from the vector GPD (cf Eq. (|5|), only lattice data exist and solely for A2${t) an d 
^2,2(i) [23]. Also for those GFF's, the comparison was quite encouraging. 

In Fig. 2, our preliminary results for the tensor GPD of a charged pion are shown 
for some values of the variable £ and t, but for < x < 1. The GPD for negative values 
of x can be obtained by exploiting the symmetry property of E^ IV (x, £, t) (see, e.g. 
Ref. [4]) under the transformation x — ¥ — x. The symmetry property follows from both 
the charge-conjugation (p — > —p and p' — > —p) and the isospin invariance and leads to 

E I £{x 1 t,t)=-E I £{-x,t,t) , E%{x,li,t) = E%{-x,t,i) (21) 

It is important to remind that for £ — > the valence component is dominant (DGLAP 
regime) while for £ — > 1 the non valence term is acting (ERBL regime). In view of 
that, it is expected a peak around x ~ 1 for £ — > 1, as discussed in Ref. Finally, 
from Eq. ( 21 1, one can understand the vanishing value of E^(0, £, t). In Ref. [TS , our 
preliminary CCQM results have been compared with the calculations obtained within 
a LF Hamiltonian Dynamics framework, where only the valence component have been 
retained |3j, obtaining a nice agreement in the relevant kinematical regions. 

From the tensor GPD, one can calculate the corresponding GFF, see Eq. |5J. We 
have started (see Ref. [18] for a preliminary presentation) the investigation by evalu- 
ating Bi } o(t) and B2 } o(t), and comparing our results with the lattice data 110, 12,23 , 
as well as with other models like the chiral quark model of Ref. [23] and the instanton 
model of Ref. |25| . The complete analysis will be presented elsewhere |20| . For a better 
comparison between our calculations of GFF's and the lattice results, we have also 
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m„ [GeV] 



Fig. 3 Preliminary B„,o(0) vs m-n, for n = 1,2. Thin solid line: CCQM Bi i0 (0,m w ). Thin 
dashed line: B* (0, m^) = 1.13 m,,-. Thick solid line: CCQM i?2.o(0, m^). Thick dashed line: 
^2 (1 (0, m T ) = 0.74 m-K- As a guide, the physical pion mass is indicated by a vertical dashed 
line. 



evaluated Bi,o(i) and E>2fi(t) both at the physical value of the pion mass, m w = 140 
MeV, and at the actual lattice value, m^* = 600 MeV, tuning the CQ mass according 
to the Chiral Perturbation Theory (see, e.g. Ref. |26|). namely 



Phys 



Physy2 



Phys 



(22) 



where m. 



Phys 



and m 



Phys 



are the constituent quark mass and the current quark mass 



at the physical pion mass, respectively. In our approach the values m^ hys — 220 MeV 



and m. 



Phys 



4 ± 1 MeV have been taken. Notice that lattice data were extrapolated 



to the physical pion mass by using a proper chiral formula |12II23) . In particular, we 
can anticipate that the behavior of Bi^(t) /Bi o(0) fully agrees with both i) the lattice 
data [12] and ii) the calculations from xQM [24] and instanton model [25], while for 
82.0 M / -B2.o(0) the comparison improves as the physical pion mass is approached. It 
is worth noting that the GFF's has been rescaled by their own value at t = 0, in order 
to get rid off the effect of the QCD evolution on the tensor GPD (cf [24] and [25]). 
Furthermore, since the tensor charges of the pion are attracting an increasing interest 
(see also Ref. [27], where -Bi,o(0) is taken as a check for the lattice evaluation of other 
observables), it is important to carefully investigate -B nj o(0), for n = 1, 2, as a function 
of a running m-n , for getting information on their chiral limit (i.e. the relevant limit for 
lattice calculations). In Fig. 3, the preliminary results of CCQM for both 61,0(0, m T ) 
and -82,0(0, jtitt) are shown; it is understood that in our analytic model, -B n ,o(0, m v ) 
is evaluated by changing the values of and m q , according to Eq. (221. It is worth 



noting that the linear behavior as a function of m^, starts for i?2,o(0, ra^) earlier than 
for Bi^O, m-jr), and in particular, for < 300 MeV one has 61,0(0, m^) ~ 1.13 m-n- 
and i?2,o(0, m w ) ~ 0.74 m^. 

To accomplish a meaningful comparison with the lattice data, one has to evolve 
the CCQM results up to the scale /1 = 2 GeV, the scale of the lattice data. According 
to [23], the evolution should affect ^2,0(0) more than 6i,o(0) (almost a factor of two 
for fiQ = 320 MeV). Though the complete investigation will be reported elsewhere [20| , 
we can anticipate here that for n = 1, the CCQM tensor charge is lower than the 
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Table 1 Mean shifts along the direction perpendicular to the u-quark transverse polarization, 
s x = {1,0}, for n = 1,2. 





CCQM - = 140 MeV 


Lattice 12 


(Mi 


0.113 fm 


0.151 ± 0.024 fm 


(by)2 


0.092 fm 


0.106 ± 0.028 fm 



lattice one presented in in Ref. |12) . that in turn is bigger than the recent lattice data 
of Ref. [27| . Furthermore, the comparison with the xQM calculation [24], in the limit 
of m-Tr — > 0, seems reasonable, since B^ M (0) ~ 1.33 and B% ? M (0) ~ 0.47 at 
the quark-model scale fio — 320 MeV. 

The calculation of the GFF's, allows one to investigate the probability density 
Pn (bj_,sj_) for a transversely-polarized u-quark (cf Eq. so that one can address 
the IPS structure of the pion. Table I shows preliminary results (still without evolution 
effects) for the average transverse shifts when the quark is polarized along the rc-axis, 
i.e. sj_ = {1,0}. The shift for a given n is 

, J dh ± by Pn(b_L,s ± ) _\_ B n , (t = 0) 
Kv)n J dh ± Pn (b ± , s± ) 2m„ A n , (i = 0) 1 ^ 

The values shown in Table I clearly demonstrate the dipole-like distortion of the 
transverse density in a direction perpendicular to the quark polarization, pointing to 
a non trivial correlation between the orbital angular momenta and the spin of the 
constituents inside a pseudoscalar hadron (see, e.g., Ref. [9] and references quoted 
therein) . 



5 Conclusions and Perspectives 

A simple, fully covariant constituent quark model has been exploited for investigating 
the phenomenology of the leading-order Generalized Parton Distributions of the pion. 
The model has been already applied to the vector GPD [HGI], an d in the present 
work some preliminary results for the tensor GPD have been illustrated. The main 
ingredients of the approach are i) an Ansatz for the Bethe-Salpeter amplitude and ii) 
the generalization of the Mandelstam formula, applied in the seminal work of Ref. [TS] 
to matrix elements of the em current operator between states of a relativistic composite 
system. 

The CCQM tensor GPD, E v j'(x,^,t), favorably compares with the one obtained 
within the Light-front Hamiltonian dynamics (see |3lll8]} and this represents an im- 
portant test of the CCQM, as well as the analysis of the tensor charges by varying 
the pion mass (see Fig. 3), relevant for the chiral extrapolation of the lattice results. 
Finally the first evaluation of the transverse shifts of the probability distribution of 
a transversely-polarized quark has been given in Table I, indicating an encouraging 
comparison with the lattice data. 

In order to construct a more and more reliable phenomenological tool, the CCQM 
study of the pion GPD will be improved by including the evolution with the energy 
scale, and by developing new Ansatzes of the Bethe-Salpeter amplitude, that allow 
one to increase the dynamical content, for instance within the Nakanishi perturbation 
integral representation of the quark-hadron vertex (see, e.g., Ref. [28] 1. 
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